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Evidence of Quasi-linear Super-Structures in the Cosmic 
Microwave Background and Galaxy Distribution 

Kaiki Taro Inoue^, Nobuyuki Sakai^ and Kenji Toinita^ 



O : ABSTRACT 



Recent measurements of hot and cold spots on the cosmic microwave back- 
ground (CMB) sky suggest a presence of super-structures on (> 100 h~^Mpc) 



^ ■ scales. We develop a new formalism to estimate the expected amplitude of 

temperature fluctuations due to the integrated Sachs-Wolfe (ISW) effect from 

O ' prominent quasi-linear structures. Applying the developed tools to the observed 

s-J ■ ISW signals from voids and clusters in catalogs of galaxies at redshifts z < 1, 

i-C ■ we find that they indeed imply a presence of quasi-linear super-structures with a 

I ■ comoving radius 100 ~ 300 h~^Mpc and a density contrast \6\ ~ O(O.l). We also 

i-{ ■ find that the observed ISW signals are at odd with the concordant A cold dark 

^ ■ matter (CDM) model that predicts Gaussian primordial perturbations at > 3 a 

level. We confirm that the mean temperature around the CMB cold spot in the 

OnJ ■ southern Galactic hemisphere filtered by a compensating top-hat filter deviates 

from a mean value at ~ 3 a level, implying that a quasi-linear supervoid or an 



lO . underdensity region surrounded by a massive wall may reside at low redshifts 

t;;-!- . z < 0.3 and the actual angular size (16° — 17°) may be larger than the apparent 

If) ■ size (4° — 10°) discussed in literature. Possible solutions are briefly discussed. 
O 



Subject headings: cosmic microwave background - cosmology - large scale struc- 
ture of the universe 



1. Introduction 

Although the A cold dark matter (CDM) models have succeeded in explaining a number 
of observations, some problems remain unresolved. For example, origins of a possible break 
of statistical isotropy in the large-angle cosmic microwave background (CMB) anisotropy 
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(Tegmark et al. 2003; Eriksen et al. 2004; Vielva et al. 2004) and a possible discrepancy 
between observed and theoretically predicted galaxy-CMB cross-correlation (Rassat et al. 
2007; Ho et al. 2008) are still not understood well. These observational results imply that 
structures on scales larger than > 100 Mpc (super-structures) in our local universe are more 
lumpy than expected (Afshordi et al. 2009). 

As the origin of the large-angle CMB anomalies, many authors have considered a pos- 
sibility that the CMB is affected by local inhomogeneities (Moffat 2005; Tomita 2005a,b; 
Cooray & Seto 2005; Rakic & Schwartz 2007). Inoue & Silk (2006, 2007) have shown that 
a particular configuration of compensated quasi-linear supervoids can explain most of the 
features of the anomalies. Subsequent theoretical analyses have shown that the CMB tem- 
perature distribution for quasi-linear structures can be skewed toward low temperature due 
to the second order integrated Sachs- Wolfe (ISW)(or Rees-Sciama) effect (Tomita & Inoue 
2008; Sakai & Inoue 2008). 

In fact, Granett et al. (2008) found a significant ISW signal at the scale of 4° — 6° at 
redshifts around z ~ 0.5 and a weak signal of negatively skewed temperature distribution for 
distinct voids and clusters at redshifts 0.4 < z < 0.75. Moreover, Francis & Peacock (2009) 
have shown that the ISW effect due to local structures at redshift z < 0.3 significantly 
affects the large-angle CMB anisotropies and that some of the CMB anomalies no longer 
persist after subtraction of the ISW contribution. These observations of galaxy-CMB cross- 
correlation may suggest an existence of anomalously large perturbations or new physics on 
scales > 100 Mpc. 

In order to evaluate the significance of the ISW signals for prominent structures, N-body 
simulations on cosmological scales seem to be suitable(Cai et al. 2010) for this purpose. 
However, the computation time is relatively long and finding physical interpretation from 
a number of numerical results is sometimes difficult. In contrast, analytical methods are 
suitable for estimating the order of statistical significance in a relatively short time, and 
physical interpretations are often simpler. 

In this paper, we evaluate the statistical significance of the ISW signals for prominent 
super-structures based on an analytic method and try to construct simple models that are 
consistent with the data. In section 2, we develop a formalism for analytically calculat- 
ing the ISW signal due to prominent non-linear super-structures based on a spherically 
symmetric homogeneous collapse model and we study the effect of non-linearity and inho- 
mogeneity of such structures. In section 3, we apply the developed method to observed data 
and calculate the statistical significance of the discrepancy between the predicted and the 
observed ISW signals. In section 4, we discuss the origin of the observed discrepancy. In 
section 5, we summarize our results and discuss some unresolved issues. In the following. 
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unless noted, we assume a concordant ACDM cosmology with {Qrn,o,^A,o,^b,o,h,crs,n) = 
(0.26,0.74,0.044,0.72,0.80,0.96), which agrees with the recent CMB and large-scale struc- 
ture data (Sanchez et al. 2009). 
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2. Cross correlation for prominent quasi-linear structures 

2.1. Thin-shell approximation 

For simplicity, in this section, we assume that super-structures are modelled by spheri- 
cally symmetric homogeneous compensated voids/clusters with an infinitesimally thin-shell. 
The background spacetime is assumed to be a fiat FRW universe with matter and a cosmo- 
logical constant A. 

Let K and ^ be the curvature and the physical radius of a void/cluster in unit of the 
Hubble radius H~^, and 6h as the Hubble parameter contrast, t denotes the cosmological 
time. We describe the angle between the normal vector of the shell and the three dimensional 
momentum of the CMB photon that leaves the shell by ip. We assume that the comoving 
radius r^ of the void in the background coordinates satisfies r^ oc t^, where /3 is a constant. 

Up to order 0{{ry/H^^)^) and 0(fi:^), the temperature anisotropy of the CMB photons 
that pass through spherical homogeneous compensated voids in the flat FRW universe can 
be written as (Inoue & Silk 2007), 

AT _ 1 
~Y' ~ 3 



r cos ij -26f, - 5^ + (3 + A6m)6H^r 
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where 6m and Qm denote the matter density contrast of the void and the matter density 
parameter, respectively. The variables ^, ip, 6m, 6h, and flm are evaluated at the time the 
CMB photon leaves the shell. It should be noted that the formula ([1]) is valid even if \6m\ 
or \6h\ is somewhat large as long as the normalized curvature n is sufficiently small. The 
formula (II]) can be also applied to spherical compensating clusters with a density contrast 
6m > corresponding to a homogeneous spherical cluster with an infinitesimally thin "wall" . 
This approximation holds only in weakly non-linear regime since the amplitude of the density 
contrast corresponding to a negative mass cannot exceed 1. We examine this approximation 
in sections 2-4 in detail. 

Because we are mainly interested in linear l^^l ^ 1 and quasi-linear l^^l = O(O.l) 
regime, we expand 6h in terms of 6m up to second order as 



JH 



nm6m{l + l/fM)/2~e6 
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where w is an equation- of- state parameter, e is a constant that describes the non-hnear effect 

and 

'5 1 11 



5 



6'3' 6 ' 



(4) 



where 2-^1 is Gauss' hypergeometric function, e can be estimated from numerical integration 
of the Friedmann equation inside the shell as we shall show later. In a similar manner, for the 
shell expansion, we assume the following relation for the wall peculiar velocity normalized 
by the background Hubble expansion, 

v = s-'P, P = -^^Sm + i^6l, (5) 

o 

where u represents a constant that describes the non-linear effect (Inoue & Silk 2007). f{^m) 
is written in terms of the scale factor a and the growth factor D as 

In quasi-linear regime, simplification e = 77 = can be verified, which will be shown in 
section 2-3 and section 2-4. 



2.2. Homogeneous collapse 

In order to describe the dynamics of local inhomogeneity, we adopt a homogeneous 
collapse model which consists of an inner FRW patch and a surrounding background fiat 
FRW spacetime (Lahav et al. 1991). The size of the inner patch is assumed to be sufficiently 
smaller than the horizon H~^ in the background spacetime. 

We assume that both the regions have only dust and a cosmological constant A. The 
time evolution of either the inner patch or the background spacetime is described by the 
Friedmann equation, 

H "m,0 ,0 , ^ ~ tot /„N 

"A,0 H o ) {') 



Hi a? ' ''" ' a2 
where a denotes the scale factor, Vtmfi, ^a.o, ^tot are the present energy density parameters of 
non-relativistic matter, a cosmological constant A, and the total energy density, respectively. 
The scale factor at present for the background spacetime is set to ao = 1. In what follows, 
we describe variables in the inner patch by putting tilde " ~" on top of the variables and we 
consider only flat FRW universes with dusts and a cosmological constant A. 

First, we calculate matter density contrast 5m of the inner patch. Initially {zi ^ 1), we 
assume that the fluctuation of the matter perturbation 5mi is so small that ai ^ ai,Hi ^ Hi. 
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Then, the Friedmann equation ([7]) yields, 

Hi' — 5^ — +^^'°"^^5r- ^^ 

In terms of physical radius of the patch R = dr, where r is the comoving radius measured 
in the background spacetime, equation ([8]) can be written as 



&J ^ -' 



(9) 



R^ 
R 



where t is the cosmological time. The matter density contrast 6m can be written as a function 
of a ratio of the present and the initial comoving radius rj = r/ri as 5m = ti~^ — 1. From 
equation ([H]), "Z] as a function of redshift z is given by solving 



drj 
dz 



= -[-Qrn,o{^ + Zi)Smi + ^mfl{v/{^ + Z)) ^ 

+ fiA,o(r7/(l + z)Y] '/'[fi^,o(l + zf + fiA,o]^'/' 



(10) 



1 + z 



Note that right-hand-side in equation flTO|) does not depend on Ri. From numerical integra- 
tion of equation flTOj) . the matter density contrast 6m as a function of redshift z is obtained 
by setting initial density contrast 6mi = Sm{zi). 

The linearly perturbed matter density contrast 6^ in the FRW background spacetime 
is given by 

^ 3 6m^H{z) [^ U + 1 

where H^{z) = fim,o(l + -z)^ + ^a,o (Heath 1977). Constant 3/5 comes from our choice 
of initial condition that the peculiar velocity inside the initial patch is zero. The relation 
between 6m and 6^ is shown in figure 1. Because the relation does not change so much even 
if one varies the cosmological parameters of the background spacetime, non-linear isolated 
homogeneous spherical patches can be solely calculated from corresponding linear perturba- 
tions (Friedmann & Piran 2001). It should be noted, however, that the relation is valid only 
if 6m ^ 6j;^cut = —0.8 because of shell-crossing (Furlanetto & Piran 2006). 

Next, we calculate the Hubble parameter contrast 6h = H/H — 1 in non-linear regime. 
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Fig. 1. — Density contrast 6^ as a function of linear density contrast 6^ for voids (left) and 
clusters (right). 
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Fig. 2. — Non-linear parameter e as a function of density contrast 6m in the EdS model. 
Plugging r] = R{1 + z)/{Ri{l + Zi)) into equation ([H]), we have 

- 5™a^,o(l + ^i)(l + 2;)V'), (12) 

where rj is given by solving equation (ITU]) . From equation (3) and (12), one can estimate the 
non-linear parameter as 



6 = Cfim,o(i + l/fiw))/2 - 6;;,\H/H - 1). 



(13) 



2.3. Effect of non-linear dynamics 

In previous section, we have seen that the non-linear density contrast 6m for a spheri- 
cally symmetric homogeneous patch can be written in terms of corresponding linear density 
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Fig. 3. — Non-linear parameter e as a function of density contrast Sm in the flat-A model 
with Qm,o = 0.26. 

contrast 6^. In order to calculate the ISW effect, we need to estimate the Hubble parameter 
contrast 6h and the peculiar velocity v of the wall. The non-linear corrections to 6h and v 
can be characterized by two parameters e and z/, respectively. 

First, we consider the effect of non-linear correction to the Hubble contrast 6h- As one 
can see in figure 2 and 3, e is always negative and the amphtude is |e| < 0.16 for l^^l < 1.0. 
This represents a slight enhancement in the expansion speed within the inner patch due to 
non-linearity. In low-density universes {flmfl < 1)? |e| is smaller than that in high-density 
universes. In the Einstein-de Sitter (EdS) universe, e depends only on 6m (figure 2). In 
contrast, in low-density universes, e depends on the amplitude of the initial epoch as well 
(figure 3). This is because the expansion speed inside the patch is suppressed when the 
energy component of the background universe is dominated by a cosmological constant A. 
We have found that the non-linear contribution to Sh is less than 10 per cent for \6m\ < 0.2 
and flm > 0.26. 

Second, we consider the effect of non-linear correction to peculiar velocity of the wall. 
In the thin-shell limit, the motion of the spherically symmetric wall can be obtained by 
numerically solving a set of ordinary differential equations using Israel's method(Israel 1966, 
Maeda & Sato 1983). If the inner region and the outer region are described by the FLRW 
spacetime, the fitting formula for the peculiar velocity of an expanding wall normalized by 
the background Hubble expansion can be written as (Maeda, Sakai, & Triay 2010) 



00.56 



+ o.isl + om\6. 



(14) 



for Qm + ^A = 1- We have confirmed that the accuracy of the fitting formula is within one 
percent for < Qm < 1 and \6m\ < 1 using numerically computed values. From equation 
flT^ . we find that the contribution of non-linear effect is less than 5 percent for \6m\ < 0.3. 
Thus an approximation 5h cc 5m ccv can be validated in the quasi-linear regime. 
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2.4. Effect of non-linearity and inhomogeneity on the ISW signal 

In literature, the thin-shell approximation has been often used to describe almost empty 
voids with (5 ~ — 1 (Maeda & Sato 1983). In quasi-linear regime, however, we also need to 
consider the effect of thickness of the wall and inhomogeneity of the matter distribution 
because quasi-linear voids are not in the asymptotic regime. Non-linearity of the wall may 
significantly affect the CMB photons that pass through it. Moreover, it seems not realistic 
to apply the thin-shell approximation to spherically symmetric clusters since the mass of the 
wall cannot be negative. 

In order to estimate the validity of the thin-shell approximation, we have compared the 
ISW signal with those obtained by using second order perturbation theory (Tomita & Inoue 
2008) and by using the Lemaitre-Tolman-Bondi (LTB) solution (Sakai & Inoue 2008), which 
yields exact results without recourse to the cosmological Newtonian approximation. We have 
assumed top-hat type matter distribution (for linear matter perturbation) for void/cluster 
for calculation using second order perturbation theory and a smooth distribution specified 
by a certain polynomial function for calculation using the LTB solution. The voids/clusters 
are assumed to be compensated so that the gravitational potential outside the cut-off radius 
(ri for the perturbative analysis and Tout ior the LTB-based analysis) is constant. For detail, 
see appendices A and B. 

As an example, we have computed temperature fluctuations AT generated from a com- 
pensated void/cluster using the three types of method. The density contrast, the comoving 
radius and the redshift of the center of a void/cluster are set to \6m\ = 0.1, r^ = 200 /i~^Mpc, 
^1 = f^out = 210 h^^Mpc, and z = 0.2, respectively. The width of the wall is assumed to 
be 1/10 of the cut-off radius. As one can see in figure 4, the three methods agree well for 
low density universes in which the linear effect is dominant. In contrast, the discrepancy 
becomes apparent for high density universes in which the non-linear effect is dominant. This 
discrepancy is partially due to a slight difference in the assumed density profile (top-hat 
type for the perturbative analysis, polynomial type for the LTB). In order to demonstrate 
the role of non-linearity, we have plotted first order (linear ISW effect) and second order 
(RS effect minus linear ISW effect) contributions to the ISW signal (figure 5). The first 
order effect makes the CMB temperature negative(positive) for a void(cluster) but the sec- 
ond order effect makes the CMB temperature negative near the center and positive near 
the boundary regardless of the sign of the density contrast. As a result, the amplitude of 
temperature fluctuation for a void(cluster) is enhanced(suppressed) in the direction near the 
center but it is suppressed (enhanced) in the direction near the boundary. These non-linear 
effects become much apparent for models with higher background density because the linear 
ISW effect becomes less effective. However, if we take into account of the thickness of the 
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Fig. 4. — Temperature fluctuation AT due to a compensated cluster (left panels) /a compen- 
sated void (right panels) centered at redshift z = 0.2 as a function of angular radius from the 
center. We used three methods for deriving the ISW signal: the thin-shell approximation 
(thick curve), the second order perturbation theory (dotted curve) and the LTB solution 
(dashed curve). We set the comoving (outer) radius of the cluster /void r„ = 200 h^^Mpc for 
the thin-shell approximation, Tout = 210 h^^Mpc for the second order and the LTB calcula- 
tions with wall width about tenth of the outer radius. The density parameter is Qmo = 0.26. 
For detail, see the text. 
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Fig. 5. — The first and the second order contributions to temperature fluctuation AT due 
to a compensating cluster (left)/void (right) centered at redshift z = 0.2 as a function of 
an angular radius from the center. We set the outer comoving radius Tout = 210 h~^Mpc, 
the inner comoving radius rin = 0.93 Tout, the density parameter Qmo = 0.6, and the density 
contrast at the center of the cluster /void \6m{z = 0.2) | = 0.1. 

wall, these non-linear effects can be less conspicuous since the amplitude of the gravitatinal 
potential becomes smaller for a flxed outer radius (flgure 6). In the A-dominated universe, a 
quasi-linear compensated void can be recognized as a cold spot surrounded by a very weak 
hot ring, whereas a quasi-linear compensated cluster can be recognized as a hot spot possibly 
with a dip at the center of it. In the EdS universe, either a compensated quasi-linear void 
or a cluster can be identifled as a cold spot surrounded by a hot ring. 



2.5. ISW effect from prominent quasi-linear structures 

In order to fully utilize information of the three dimensional distribution, we consider a 
temperature anisotropy AT/T obtained from stacked images on the CMB sky that corre- 
sponds to most prominent voids/clusters in a galaxy catalog. First, we flx an angular radius 
6out of a circular region on the CMB sky that will be used in the stacking analysis. Then, the 
corresponding smoothing scale r^ in comoving coordinate for the corresponding fluctuation 
at z is Tg = {1 + z)DA{z)9outi where Da{z) is the angular diameter distance to the galaxy. 
The corresponding initial smoothing scale is r\ = rj~^rs. 

We assume that the probability distribution function (PDF) of linear density contrast 
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Fig. 6. — The effect of thickness of the wall for Q^no = 0.6 ri = 210/i~^Mpc, and 6m{z = 
0.2) = 0.2 (left) and 6m{z = 0.2) = —0.2 (right) for top-hat type density perturbations (see 
appendix A). 
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Fig. 7. — Effect of non-linearity: the PDF of Sm (dashed curve) and that of S^ (full curve) 
for a smoothing scale r^ = 50/i~^Mpc at z = 0.5. 



(5^ at redshift z is given by a Gaussian distribution function, 
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where o^{r\^z^ is the variance of the linearly extrapolated density contrast at redshift z 
smoothed by a spherically symmetric top-hat type window function with an initial comoving 
radius r\. Note that cr(r*,z) depends on cosmological parameters such as r2m,o,cr8 and the 
spectrum index n. Then, the PDF of non-linear density contrast of the inner patch 6^ is 
given by 

P^^(5^; a{T\, z)) = aP\5t{5^)- a{rl, z))-^, (16) 



d6„ 



where a is a constant that normalizes the PDF. 



As shown in figure 7, the PDF of 6m is positively skewed in comparison with the PDF of 
6m because of non-linearity. For a sample region at redshift z with a total comoving volume 
V, the total number of voids or clusters with a radius r^ is approximately Nt ~ 3V/{A7irl). 
In what follows we assume that the number of prominent voids/clusters {Ny/Nc) determines 
the corresponding threshold of density contrast 6m,th{z), which is given by 



and 



NjNt 



NjNt 



P^'^{6m;cT{rl,z))d6„ 



(17) 
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respectively. From equations ([T]), flTTl) . and flTSl) the mean temperature fluctuation within an 
angular radius Oout for a stacked A^^, or N^. images corresponding to prominent quasi-linear 
voids/clusters at redshift ~ z can be approximately written as 

^ ^ 7r9iJP^^iSm;a)dSm ' ^ ^ 

where < ^ < Oout, W{6; 6in) is a compensating window function that satisfies 

/ 2TTew{e; ein)de = - 2TTew{e; ein)de (20) 

Jo J din 



and 

ip{9) = 9 + arcsin 

where h = rs{l + z)/Da{z). 



h ^ cos 9 sin 9 — v sin^6' — /i~^ sin'^6' 



(21) 



3. Application to observations 

3.1. SDSS-WMAP cross correlation 

A cross correlation analysis using a stacked image built by averaging the CMB sur- 
rounding distinct voids/clusters has been done by Granett et al. (2008). They have used 1.1 
million Luminous Red Galaxies (LRGs) from the SDSS catalog covering 7500 square degrees. 
The range of redshifts of the LRGs is 0.4 < z < 0.75, with a median of ~ 0.5. The total 
volume is ~ 5/i~^Gpc^. They used so-called the ZOBOV (ZOnes Bordering On Voidness; 
Neyrinck 2008) algorithm to find supervoids and superclusters in the LRG catalog and made 
a stacked image from an inversely variance weighted WMAP 5-year (Q,V, and W) map. In 
order to reduce contribution from CMB fluctuations on scales larger than the objects, they 
used a top-hat type compensating filter 

where 9 out = cos"^ (2cos6'i„ — 1). 

First, using the developed tools based on thin-shell approximation and homogeneous 
collapse model in section 2, we estimate the expected amplitude of the ISW signal for promi- 
nent structures in a concordant ACDM model with Gaussian primordial fluctuations and 
compare with the observed values obtained from the SDSS-LRG catalog. The number of 
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most distinct voids or clusters A^ and the cut-off radius din are chosen as free parameters. 
At redshift z = 0.5 , the mean density contrast filtered by a top-hat type function with 
radius r = 130 h~^Mpc corresponding to 9out ~ 5.6° is just (6^) = 0.046 and the background 
density parameter is Qm{z = 0.5) = 0.54. Because the influence of non-linear ISW effect 
is weaker than that of the linear ISW effect in this setting, we expect that details of non- 
linear calculations will not much affect the result. In what follows, we use an approximation 
6h cc 6m cc V, where 6m is determined from the homogeneous collapse model in section 2. 

As shown in table 1, it turned out that the expected values of the ISW(Rees-Sciama) 
signal are typically of the order of O(10~^)K. As expected, the amplitude gets larger as 
the number of stacked image decreases, and the amplitude for voids systematically becomes 
larger than those for clusters by 5 — 10 percent (Tomita & Inoue 2008; Sakai & Inoue 2008). 
On the other hand, the order of the observed amplitudes are extremely large as O(10^^)K. 
It turns out that the discrepancy remains at 3 — 4(t level for A^ = 30 and A^ = 50. 

Second, we reconstruct the mean density proflle from the observed ISW signal for the 
SDSS-LRG catalog using our LTB model. From figure 7, one can notice a hot ring around a 
cold spot for the stacked image of voids and a dip at the center of a hot spot for the stacked 
image of clusters. Although the amplitude of the hot-ring cannot be reproduced well, the 
observed dip at the center of the hot spot can be qualitatively reproduced in our LTB models. 
We have found that the dip at the center of a compensated cluster can be generated only 
if the linear ISW effect balances the non-linear ISW effect in a limited parameter region. 
Thus the observed features in stacked images strongly imply that the corresponding super- 
structures are not linear but quasi-linear or non-linear objects. The density fluctuations 
which are necessary to produce the observed ISW signals are found to be tremendously 
large. In flgure 8, we plot the ISW signal from a compensated cluster with 6m ^ 7 o- and 
that from a compensated void with \6m\ ~ lOo" at 2; = 10 (see the radial density proflles at 
flgure Al in appendix A). Even for these very rare objects, the amplitudes of ISW signals in 
our LTB models are much smaller than the observed ones. In fact, the mean temperatures 
for a compensating fllter 6in = 4° are 3.6/iK(1.4 a) and — 3.1/iK(2.6 a) for the cluster and the 
void, respectively. On the other hand, the probability of generating these fluctuations is as 
extremely small as 10"^^ in standard inflationary models that predict primordial gaussianity. 

Thus, the observed large ISW signals for the stacked image strongly suggest a presence 
of super-structures on scales 0(100 /i~^Mpc) with anomalously large density contrast 0(0. 1) 
which can not be produced in the concordant LCDM model. 
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Fig. 8. — Observed mean CMB temperature radial profiles for a stack image of 50 
clusters (left) and that of 50 voids (right) (thick full curves) in the SDSS-LRG cata- 
log (Granett et al 2008) and the corresponding theoretical radial profiles for compen- 
sated spherical clusters (left) /void (right) based on the LTB solution (dashed curves) and 
the linear calculation (thin full curves) with a top-hat type compensated density distri- 
bution corresponding to the LTB solutions. The LTB parameters for the set of clus- 



ters are {r out, Tin, 6m, z^) = (0.055ifo 



, 0.022i7o"', 



1.0,0.5) and those for the set of voids 



are (0.050ifQ" ,0.015ifQ" ,—0.55,0.5). The effective radii of the inner patch at a redshift 
2; = 10 and the mean filtered temperature with 6'j„ = 4° are {0.029Hq^, 3.6/iK) (cluster) and 
{0.023Hq^, — 3.1/iK) (void), respectively. Evolution of the density profile in the LTB models 
are shown in appendix B. la error bars are obtained from 1000 Monte Carlo simulations on 
the WMAP7 Q+V+W map smoothed at 1° scale (see section 3.3). 



3.2. 2MASS-WMAP cross correlation 



Francis & Peacock (2009) estimated the local density field in redshift shells using pho- 
tometric redshifts for the 2MASS galaxy catalog. They reconstructed the CMB anisotropies 
due to the ISW effect from the local density field 6m- They approximated the bias in each 
redshift shell by a linear bias relation 6g = b6m and assumed that the bias is independent 
of scale and redshift in each shell. In order to obtain the bias parameter b, a maximum 
likelihood analysis of the galaxy catalog was performed. 

There are two prominent spots in the reconstructed CMB anisotropy. One is a hot 
spot due to a supercluster around the Shapley concentration at redshifts 0.1 < z < 0.2. 
Another one is a cold spot due to a supervoid at redshifts 0.2 < z < 0.3 in the direction to 
(/, b) ~ (0, —30°). The angular radii of both structures are 9 = 20° — 30°. The temperatures 
near the center of both structures are ~ 20 fiK. The position of the supervoid is very close 
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Fig. 9. — Left: Theoretically modeled temperature profiles for an observed (full curves) and 
an expected (dotted curves) supercluster in the 2MASS catalog. Right: modeled temperature 
profiles for an observed (full curves) and an expected supervoid (dotted curves). These 
superstructures reside at redshifts z = 0.1 — 0.3. We have plotted two possible profiles for 
each structure since there is an ambiguity in the angular size due to errors {Az ~ 0.03) in 
photometric redshifts in the 2MASS galaxy catalog (Francis & Peacock 2009). 

to the one predicted in Inoue & Silk (2007), (/, 6) ~ (330°, -30°). 

Based on developed methods in section 2, we have estimated the expected density 
contrast 6m and the corresponding temperature profile (figure 9) due to a most prominent 
object in the shell. We have assumed the same cosmological parameters and primordial 
gaussianity as those discussed in section 3. In order to compute the temperature profile, 
we have used a homogeneous thin-shell model. As shown in table 2, the observed density 
contrasts are larger by 4-7 times the expected values. If the comoving radius of the structure 
is ~ 200 h~^Mpc, the absolute value of the density contrast should be \6m\ = 0.2 — 0.3, which 
implies a presence of anomalous quasi-linear super-structures. Our result is consistent with 
the power spectrum analysis in Francis & Peacock (2009) where a noticeable excess of the 
observed power at low multipoles 2 < / < 4 was reported. 



3.3. The CMB cold spot 

The most striking CMB anomaly is the presence of an apparent cold spot in the Wilkin- 
son Microwave Anisotropy Probe (WMAP) data in the Galactic southern hemisphere (Vielva 
et al. 2004; Cruz et al. 2005) (see figure 10). The cold spot has a less than 2 per cent prob- 
ability of being generated as random gaussian fluctuations (Cruz et al. 2007a), if one uses 
spherical mexican-hat type wavelets as filter functions (see also Zhang & Huterer 2010). 
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Table 1. Expected and observed amplitude of mean temperature for a compensating filter 



e^n = 4° 



A^ cluster (fiK) void (fiK) average (/i-ft') 



1 


0.98 






-1.2 




1.07 


5 


0.82 






-0.94 




0.88 


10 


0.73 






-0.83 




0.78 


30 


0.57 






-0.64 




0.61 (11.1 ±2.8)^^ 


50 


0.48 


(7. 


9 ±3.1)^ 


-0.53 (- 


•11.3 ±3.1)^ 


0.51 (9.6 ±2.2)^ 


70 


0.42 






-0.46 




0.42 (5.4 ±1.9)^ 



^Taken from Granett et al. (2008). 



Table 2. Expected and observed density contrast for super-structures in 2MASS galaxy 

catalog 



radius** 


expected 


observed 


radius** 


expected 


observed 


230 
150 


0.037 
0.094 


0.20 
0.69 


370 

250 


-0.013 
-0.037 


-0.049 
-0.15 



*The unit of the radii is h ^ Mpc. 
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Fig. 10.— Left: the WMAP7 ILC temperature map (40° x 40°) smoothed at 1° scale. Right: 
the averaged radial profile of the ILC map as a function of inclination angle 6 from the center 
of the cold spot (/, b) = (207.8°, —56.3°). A peak at 6' ~ 15° corresponds to a hot ring. 

Assuming that it is not a statistical artifact, a variety of theoretical explanations have 
been proposed, such as galactic foreground (Cruz et al. 2006), texture (Cruz et al. 2007b), 
and Sunyaev-Zeldovich (SZ) effect. However, these models failed to explain other large-angle 
anomalies by the same mechanism. 

Inoue & Silk (2006,2007) proposed that the cold spot may be produced by a supervoid at 
2; < 1 in the line-of-sight due to the ISW effect and have shown that another pair of supervoids 
that are tangential to the Shapley concentration can explain the alignment between the 
quadrupole and the octopole in the CMB. Subsequently, Rudnick et al. (2007) found a 
depression in source counts in the NRAO VLA Sky Survey(NVSS) in the direction to the 
cold spot, although the statistical significance has been questioned (Smith & Huterer 2010). 
Recent optical observations (Granett et al. 2009, Bremer et al. 2010), however, revealed 
that any noticeable supervoids at 0.35 < 2; < 1 in the line-of-sight are ruled out. These 
observations suggest that the angular size of the supervoid may be larger or smaller than 
expected and that it resides at low redshifts z < 0.35 or at high redshifts z > 1. 

In order to test such a possibility, we have calculated averaged temperature around the 
cold spot (see figure 11) using a spherical top-hat compensating filter Wth{0; Om). 

Interestingly, we have discovered two peaks in the plot of the filtered mean temperature 
around the cold spot as a function of inner radius of the filter (figure 12). The inner and the 
outer peaks are observed at ^i„ = 4°-5°(^o„i = 6° -7°) and Oin = 12° -13° {eout = 16° -17°). 
The outer peak corresponds to a hot ring, which is visible by eyes (see figure 10). 



In order to estimate the statistical significance of the peaks, we have used a WMAP 
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Fig. 11. — Left: Icr of the filterd mean temperature as a function of an inner radius Bin for 
the Q+V+W map (dashed curve) and that derived from Q's (thin curve) and from pseudo- 
C;'s (thick curve). Right: la for the ILC map (dashed curve) and corresponding theoretical 
values (thin curve and thick curve) as shown in the left figure. 



7-year internal linear combination (ILC) map smoothed at 1° scale with a Galactic skycut 
|6| < 20° and a combination of the Q, V, and W band frequency maps smoothed at 1° scale 
averaged with weights inversely proportional to the noise variances with a "standard" Galac- 
tic skycut made by the WMAP team. In order to reduce possible residual contamination 
from the Galactic foreground, we further cut a region |6| < 35° for the Q-l-V-l-W map. In or- 
der to estimate the errors, firstly, we generated 1000 random positions on the ILC(|6| > 20°) 
and the Q-fV+W(|6| > 35°) maps, and then computed variances cr^ for the filtered mean 
temperature. Second, we have calculated expected a^ for the filtered temperature using the 
angular power spectrum Ci for the WMAP 7-year data obtained by the WMAP team (see 
appendix A). Note that we have computed pseudo-C/'s from the Q' for each skymap. As 
shown in figure 11, the observed standard deviations are Icr = 19 ~ 20 /iK for 6'j„ = 4° and 
lo" = 14 ~ 16 /iK for 6in = 14°. Our result for 6'j„ = 4° is roughly consistent with the values 
for stacked images in Granett et al. (2008) assuming no correlation between voids/clusters 
in a particular configuration. In the Q+V+W map, a slight suppression in a is observed 
at 6in > 8°. Theoretically calculated values are found to be systematically lager than the 
observed values by 4 — 16 per cent for 4° < 6in < 14°. These discrepancies represent an 
uncertainty due to the Galactic foreground emission. 

As one can see in figure 12, a deviation corresponding to the inner peak is roughly 4 a 
and that corresponding to the outer peak is 4 ~ 4.5 a. Assuming that the filtered mean 
temperature obeys Gaussian statistics, the statistical significances are P{> 4cr) = 6 x 10"'^ 
and P{> 4.5 a) = 7 x 10"^ per cent. The total solid angle of the ILC map (|6| > 20°) is 8.27 
sr and that of the Q+V+W map (|6| > 35°)is 5.36 sr. The total number of the independent 
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Fig. 12. — Left: mean temperature around the center of the cold spot for a compensating 
top-hat filter (full curve) and for an uncompensated top-hat filter without a ring (dashed 
curve) as a function of the inner radius 6^. Right: mean temperature around the center of 
the cold spot divided by the standard deviation for the ILC map (dashed curve) and the 
Q-l-V+W map (full curve) as a function of Bin 



patches is roughly given by a ratio between the solid angle of the map and the area of the 
spherical patch with angular radius 9out- Therefore, for 9out = 7°, we have ~ 180 samples for 
the ILC and ~ 110 samples for Q+V+W map, yielding (110-180) x P{> 4 a) = 0.7-1 per 
cent. In a similar manner, for 6out = 17°, one can easily show that the statistical significance 
is ~ 0.01 — 0.2 per cent corresponding to ~ 3 a if the likelihood function is a Gaussian one. 

Thus, the cold spot surrounded by a hot ring at scale ~ 17° is more pecuhar than the 
cold spot at scale ~ 7°. Therefore, the real size of the supervoid is expected to be larger than 
the apparent size of the cold region. Because 1 a deviation (corresponding to 15 — 20/iK) 
due to a supervoid is enough to make the signal non-Gaussian, it is reasonable to assume 
that the contribution from a supervoid is less important than that from other effects due to 
acoustic oscillation or Doppler shift at the last scattering surface. For instance, a supervoid 
with a density contrast 6m = —0.3 with a comoving radius r = 200 h'^Mpc at a redshift 
z = 0.2 corresponding to an angular radius 9 = 20° would yield a temperature decrease 
AT ~ 20/iK in the direction to the center. Moreover, if the supervoid is not compensating, 
a wall surrounding the supervoid could generate the observed hot ring. It may consist of just 
an ordinary underdense region surrounded by massive superclusters. Further observational 
study is necessary for checking the validity of the "local supervoid with a massive wall" 
scenario. 
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4. Possible solutions 

Why the observed ISW signals for prominent structures are so large? 

One possible systematic effect may come from a deviation from spherically symmetric 
density profile that we have not considered. Indeed, gravitational instability causes pan- 
cake or needle like structures in high density regions. However, as we have seen, the order 
of the density contrast of relevant prominent super-structures is 5m = O[0.1]. Therefore, 
we expect that the effect of anisotropic collapse plays just a minor role. Moreover, in the 
case of supervoids, a deviation from spherical symmetry is suppressed as the void expands 
in comoving coordinates. Thus, it is difficult to attribute the major cause to the deviation 
from spherical symmetry. 

Another possible systematic effect is our neglect of fluctuations on larger scales. For 
instance, we may have observed just a tip of fluctuations whose real scale extends to r > 
1000 /;,~^Mpc. Indeed, the amplitude of the ISW effect is roughly proportional to the scale 
of fluctuations, i.e., AT/T oc r for r > 100/i~^Mpc (Inoue & Silk 2007). Therefore, the 
observed large amplitude of ISW signal can be naturally explained. However, the angular 
sizes of the observed hot and cold spots from the stacked images are just 4° — 6° at 2; ~ 0.5 
corresponding to r = 100 — 140 /i^^Mpc. It is difficult to explain why the angular sizes are 
so small since contributions from the ordinary Sachs- Wolfe effect and the early ISW effect 
generated near the last scattering surface are significantly suppressed by stacking a number 
of images. 

Then what are the possible mechanisms that can explain the anomalously large ISW 
signals? 

One intriguing possibility is that the primordial fluctuations are non-Gaussian. Our re- 
sults suggest that the number of both supervoids and superclusters is significantly enhanced 
in comparison with the standard Gaussian predictions. Therefore, the effect of deviation 
from Gaussianity may appear in the statistics of 4-point correlations in real space or trispec- 
trum in harmonic space. It can be also measured by the Minkowski functionals that contain 
information of 4-point or higher order correlations. At the last scattering surface, the co- 
moving scale of 300 /i~^Mpc corresponds to angular scale ~ 2°. If the background universe is 
homogeneous, such a non-Gaussian feature must appear at the CMB anisotropy at multiple 
/ ~ 100 corresponding to angular scale ~ 2° as well. However, so far no such a noticeable de- 
viation from Gaussianity in the CMB anisotropies has been observed (Vielva & Sanz 2010). 
Therefore, it is difficult to explain the observed signals by a simple non-Gaussian scenario 
unless one gives up the cosmological Copernican principle (Tomita 2001). 

Another possibility is a certain feature on the power spectrum of primordial fiuctuation 
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(Ichiki et al. 2009). Spike-like features in the primordial power spectrum appear in some 
inflationary scenarios that produce primordial black holes (Ivanov et al. 1994, Juan et al. 
1996, Yokoyama 1998). Although there is no natural reason to have a feature only on the 
scale of super-structures (/ ~ 100), observational constraints are not stringent since one 
needs to increase the number of samples if one abandons an assumption of the smoothness 
of the primordial power. 

Some cosmological models containing time evolving dark energy/quintessence or those 
based on scalar-tensor gravity predict an enhancement in the ISW effect due to an enhance- 
ment in acceleration of the cosmological expansion or non-trivial time evolution of dark 
energy or scalar field that may couple to matter or metric (Amendola 2001, Nagata et al. 
2003). This may help to explain the anomalously large ISW signal. However, at the same 
time, we need to suppress the ISW contribution on large angular scales since the observed 
angular power of the CMB anisotropy at very large angular scales / ~ 2 is relatively low. 
Models based on some alternative gravity might be helpful for realizing these observational 
features (Afshordi et al. 2009). 



5. Conclusion 

In this paper, we have shown that recent observations imply a presence of quasi-linear 
super-structures with a comoving radius 100 — 300 /i^^Mpc at redshifts z < 1. Observa- 
tions are at odd with the concordant ACDM cosmology that predicts Gaussian primordial 
perturbations at > 3 o" level. 

First, we have developed a formalism to estimate the amplitude of the ISW signal 
for prominent structures based on thin-shell approximation and the homogeneous collapse 
model. From comparison with other calculations based on perturbation theory and the LTB 
solution, we have found that our simple model works well for estimating the ISW signal for 
quasi-linear superstructures in A-dominated universes. 

Secondly, we have applied our developed tools to observations of the ISW signals using 
the SDSS-LRG catalog, the 2MASS catalog, and the cold spot in the Galactic southern 
hemisphere in the WMAP data. The ISW signals from stacked images for the SDSS-LRG 
catalog is inconsistent with the predicted values in the concordant ACDM model at more 
than 3 a. The radial profiles of the stacked image show a hot-ring around a cold spot for 
voids and a dip at the center of a hot-spot without a cold-ring for clusters. These non-linear 
features are also reproduced by our models using the LTB solutions although the agreement 
is not perfect. The asymmetrical features suggest that the observed super-structures are in 
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quasi-linear regime rather than hnear regime. The amphtudes of the ISW signals obtained 
from the 2MASS catalog at redshifts 0.1 < z < 0.3 are found to be several times larger 
than expected values. We have confirmed that the mean temperature around the cold spot 
filtered by a compensating top-hat filter with angular scale 6out = 16° — 17° deviates from 
the mean value at roughly 3 a level suggesting a presence of a hot-ring around the cold spot. 
Note that our finding is consistent with the previous result that the cold spot itself is not 
unusual but the hot-ring plus the inner cold region is found to be unusual (Zhang & Huterer 
2009). This implies that a supervoid may reside at low redshift z < 0.3 and the angular size 
may be larger (= 16° — 17°) than considered in literature (Masina & Notari 2009). 

Finally, we have discussed possible causes of the discrepancy between the theory and ob- 
servation, namely, observational systematics, primordial non-Gaussianity, features in power 
spectrum, dark energy or alternative gravity. 

We have not considered effects of non-spherical collapse which are important for im- 
proving estimation of the mass function of non-linear objects and effects of uncompensated 
mass distribution for super-structures. The extension of our analysis to more elaborate ones 
incorporating these effects would be helpful for realizing detailed comparison between the 
theory and the observation. 

Future surveys of the CMB, galaxy distribution, weak lensing and theoretical studies on 
dark energy/alternative gravity and inhomogeneous cosmology will certainly yield fruitful 
results for solving the puzzle. 
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A. First-order and second-order ISW effects 

In what follows, we derive analytic formulae for computing temperature fluctuations due 
to the ISW effect for spherically symmetric compensated top-hat type density perturbations 
using first-order and second-order perturbation theory (Tomita & Inoue 2008, abbreviated 
as TI hereafter). The relation between density perturbations of the growing mode and the 
potential function F of spatial variables are given in equations (3.6) and (3.11) of TI. The 
top-hat type density profile is parametrized in terms of two constants h and c, representing 
first order density contrasts at the center and at the wall at the present time (figure Al). 

The first-order density contrast 5^ at a conformal time t] when the scale factor is equal 
to a{ri) can be written in terms of the background matter density ps and the growth function 
P{ri) corresponding to the growing mode of density perturbation as 

5lk = ^P^'^/PB = -^,[ia'/a)P' - l](c, -6), (Al) 

for (r < ro, ro < r < ri), respectively, where a prime denotes a partial derivative with respect 
to conformal time f]. 

The second-order density contrast is expressed as 
Si ^ 6p^'^/pn = ^{(Cl-^C2)c^ [2ic + bnro/rr + b']C, + hc + bnro/rr-b']C2} (A2) 

for (r < ro. To < r < ri), respectively, where Ci and (2 are given in equation (2.19) of TI. Here 
we have omitted the terms that are negligible if ri <^ H~^ because we assume that typical 
size of super-structures is O{100)h^^Mpc. Neglecting the terms higher than second-order, 
the total density contrast Sp/pB can be written as 

PB 1 1 rL , /jLn2 



1 - C + {& - €■ (A3) 



PB + 5p 1 + 5p/pB 

For a central value of total density contrast, {Sp/p)c, we have the relation 

a{z)c^ + l3{z)c--f{z) = 0, (A4) 

where z is the redshift, a{z) = |(Ci - |C2)/(PBa^) - Wiz)]"^, f^{z) = {^p' - l)/{pBa^) and 
^{z)^l-l/[l + {5p/p)l 

In the text we consider models of supervoids and superclusters with a given set of tq, ti 
and {6p/p)c{zc), where {6p/p)c{zc) is {6p/p)c at the epoch of redshift Zc- From this set we 
obtain c and b, solving the above equation as 

c = -[(3{z,) + ^(3{z,y + 4a{z,)-f{z,)]/{2a{z,)), (A5) 
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and b is related to c as b/c = l/[(ri/ro)^ — 1] for compensated super-structures. 

The first-order and second-order temperature fiuctuations AT'^'^^T and AT^'^'> /T are 
defined by equations (4.2) and (4.4) of TI. Their expressions for a hght path passing the 
center of spherical voids and clusters are given in equations (5.11) and (5.13) of TI. For the 
other light paths, the first-order temperature fiuctuation is derived from equations (5.8) and 
(5.9) with equation (C6) of TI and expressed as 

AT«/T = -2{- + [- - 3{-f]P'}c{rofJ{r/ro), (A6) 

a a a 

where J{r/ro) is given in equation (C7) of TI for r < tq. For ri > r > tq, we have 

J{u) = -]:2ui^[h + I3- 3ui^l2]/{ui' - 1), (AT) 

o 

where u = r/ro, Ui = ri/rg and 

h = f \u^- uT^'^du = ln[(Mi + ^Jul-u-')lul 

J u 

/•Ml 



/3 = / {u'-uY'/\^du=^{u^ + 2u^){m'-uY'- (A8) 

Ju 3 

The second-order temperature fiuctuation can be derived from equations (2. 17), (2. 18), (4. 4) 
and (4.5) of TI and expressed as 

AT^^)/T = -[C[ / dXiFrf + C2 X 100 / d\%], (A9) 

Jo ' Jo 

where J^ d\{F^rY ^"^^ /o d\<^Q for r < r^ are given in equations (C3) and (C4) of TI and 
the expression of (,[ and (^2 is shown in equations (4.6) and (4.7) of TI. For ri > r > Tq, we 
have 

/•oo 1 

Jo y 

r°° 1 

100/ <l>orfA = -c\r^f[ui''h-2h + <dui^h-SuiHi]/{u{'-l)\ (AlO) 

Jo 4 

where Jj (? = 1 — 3) are given above and /4 is 

h= I \u'-u^)-^'^u-^du = -l-{-tan-^Uui^-uyu] + \/ u^' - v? I u^"} . (All) 
Ju 2u^ u 

When we compare the temperature fiuctuations in the perturbative model (in appendix 
A) and those in the full non-linear model (in appendix B), we should notice the difference 
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Fig. 13. — (Appendix A,B) Left: matter density contrast (at linear order) for a top-hat 
type spherical void used in our perturbative analysis. Right: density profiles 6m in our LTB 
models as a function of physical radius in unit of present Hubble radius for a compensated 
cluster and void. The initial condition is set at a redshift z = 10. We assume a concordant 
FRW cosmology with {flm,o, ^a,o) = (0.26, 0.74) as the background spacetime. 



of their density profiles, i.e. the top-hat profile (in appendix A) and the Sakai-Inoue (SI) 
profile (in appendix B). For our comparison in this paper, we simulate the top-hat profile 
to the SI profile by equating the outer boundaries and their zero points as follows. Here we 
represent the SI profile using the radial coordinate r defined in the perturbative model. In 
the top-hat profile, the radii in the outer boundary and the zero point are r = ri and tq, 
respectively, and in the SI profile the radius in the outer boundary is r = Vout and the zero 
point is r = (rj„ + rc)/2 approximately, in which Tc = {rout + rin)/^. If we equate these outer 
boundaries and zero points, we obtain 

ro = (r»„ + r,)/2 = (S^^ + ro„t)/4. (A12) 

Then for relative widths Wth = 1— tq/ti and wsi = ^—rin/rout, we have a relation wth = jWsi- 
In the text we show the temperature fluctuations in both models with parameters which 
satisfy this relation. 
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B. Method of computing Rees-Sciama effects using the LTB models 

Any spherically symmetric spacetime which includes dust of energy density p{t, r) and 
a cosmological constant A can be described by the LTB solution, 

ds' = -de + ^-^dr^ + R\t, r){de' + sin^ Odif^) (Bl) 

l + /(r) 

which satisfies 

R'-'-^ + ^R' + m (B2) 

_ m'jr) 

where ' = d/dr and ' = d/dt. 

Our model is composed of three regions: an outer fiat FRW spacetime, an inner nega- 
tively/positively curved FRW spacetime, and an intermediate shell region described by the 
LTB metric. At the initial time t = ti, which we choose as Zi = 10, we define the radial 
coordinate as R{ti,r) = r, and we assume (figure Al) 

1 + 5m for r < Tin 

m{r) = ^7rR'pmW{r), W{r) = { 1 + ^(8 - 15X + lOX^ - SX^) for nn<r<rout 

1 for r > rout 

(B4) 
where 

^Y* ly v* I ly 

X' ' c ' out "I'm /T-i r \ 

= 77-5 fc = , W = Tout - Tin. (B5) 

w/2 2 

Initial velocity field, v = R — HR, is given by the linear perturbation theory (TI). Then /(r) 
is determined by equation (IB2p . Our model parameters are ^2^,0, ^out, w, the redshift of the 
center of a void/cluster, Zc, and 6m{zc). 

The wave 4-vector k^ of a photon satisfies the null geodesic equations, 

F^— , k^k, = o, ^ + r(;,rr = o (B6) 

For null trajectories on the 6 = n/2 plane, the geodesic equations flB6l) with the metric fIBll) 
reduce to 

{k'f = X^ + R\k'^f, X= /^' r, (B7) 

V 1 + / 



R^k'P = constant, (B8) 



-29- 
dt fc*' dt ¥' ^ ' 

We use the null condition flB7p not only to set up initial data but also to check numerical 
precision after time-integration. 



To integrate the geodesic equations ( ]B9|) - (1B11|) together with the field equations ( ]B2p 



and (IB 131) numerically, we discretize r into A^ elements, 

ri = iAr, i = l,...,N, Ar = ^. (B12) 

and any field variable $(t, r) into $i(t) = $(t, rj). Evolution of Ri(t) is determined by 
flB2l) . but we also need data of R^it) and R'-{t). Because finite difference approximation, 
R'iit) ^ (-Rj+i — -Rj_i)/(2Ar), include errors of O(Ar^), we numerically integrate with time. 



-^ ^R + /' + -ARK 

2R\ R R^ ^3 J 



which is given by differentiating (IB2p with respect to r. Furthermore, to vanish R"(t,r) in 
the geodesic equations, we have introduced an auxiliary variable X. To prepare geometrical 
values between grid points r^ and rj+i, we adopt cubic interpolation: at each time t* any 
variable $(t*, r) in Tj < r < Tj+i is determined by 





Hu 


, r) = ax^ 


+ bx^ + ex 


+ d, 


X 


Ar 
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+ 9($, + <l>i+i) 


-$i+2 



24Ar ' 16 



(B14) 



We also have to consider null geodesies from an observer to the void/cluster. Suppose 
that the observer is at the origin and the center of the void/cluster is located at x = Xc on 
the y axis. Then, without loss of generality, on the x-y plane we can analyze light rays which 
reach the observer. Some position on the outer shell and the four momentum of the light 
there in the observer-centered coordinate are given by 

X = Xc + Tout cos (fi, y = routsm(p, (B15) 
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where E is the photon energy and a is the angle between the hght ray and the x-axis. 
Defining li^z) as a co moving length from the observer to the photon, we can write the light 

path as 

r dz 
x = l(z) cos a, y = l(z)sma, Hz)= / — . (BIT) 

Jq aoH{z) 

The solution of (JBlSp and (JBlTp gives 

I = Xccosa — \/r^^ — x'^sin'^ a, siny9 = sin a, (B18) 

^ rout 



and the null vector in the void/cluster-centered spherical coordinate, 

at the time when the photon leaves the shell, -Zieave- 

Our computing algorithm is summarized as follows: 

(i) Suppose f2m,o; ^m{zc), Tout, w and rcdshift of the ccutcr of a void/clustcr, 2;c. For each 
angle, a, Eqs. OBlSp and (1B19P give "initial" conditions of the null geodesic at z = ^leave 
in the void/cluster-centered spherical coordinate. 



(ii) Solve the field equations in the LTB spacetime, flB2p and fIBlSp . from z = Zi to z^ 



(iii) Solve the null geodesic equations (IB7P -f pTT]) together with the the field equations f lB2 



and ( 1B13P backward from z = zieave to the time when the photon enters the shell. 



C. Temperature Variance for top-hat compensating filter 

In what follows, we derive analytic formulae for computing variance of temperature 
fiuctuations on a sky for a circular top-hat compensating filter Wth{0', Oin)- We assume that 
an ensemble of the CMB fiuctuations can be regarded as an isotropic random field on unit 
sphere S"^. Let AT(6', 0) be a temperature fiuctuation at spherical coordinates {0,(f)). Then 
filterd temperature fiuctuation centered at the "north" pole {6 = 0) can be written as 

/•2tt rdout 

ATf = A'' d<p desmeWthie;e,n)AT{e,(i)), (ci) 

Jo Jo 
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where A = 27r(l — cos^j„). Plugging AT expanded in spherical harmonics Yim, 

AT = Y,ciimYi^ (C2) 

Im 

into equation (Al), we have 

AT/ = A-i5^a,oG'z, (C3) 

I 

where 



v/vr(2/ + 1) 



'2{~XinPl{Xin) + Pl~l{Xin) j + X outPl\-C out) ~ Pl-li^out) 



/ + 1 
•^in COS (7j,2? -^out '^•^in -L. V^^J 

Note that we have used a formula for the Legendre function P;, 

dPi{x) l{l + l) [^ 



dx 1 — x^ 



Pi{x)dx (C5) 



in deriving equation A4. Because AT is assumed to be isotropic on 5^, the variance of ATj 
can be written as a function of the angular power spectrum Ci as 

a) = A--Y.C,Gl (C6) 

If the CMB sky is smoothed by a Gaussian beam with the FWHM Og., then the variance 



IS 



a) = A-^Y.^,B,Gl (C7) 



I 



where Bi = exp[-a^J{l + l)/2], and a, = (81n2)-i/2^,. 

In the absence of complete sky coverage, we cannot directly observe Ci. We can only 
compute estimated expansion coefficients for the observed region R in the sky (Bunn 1995), 



km = Ni„, f ATYir^ sin 9 d9d<t>, {CI 

JR 



where Ni^ is a factor chosen to normalize him appropriately. If R is azimuthally symmetric, 
one possible prescription is to set (Peebles, 1973) 

Nlra = W-^\ (C9) 

where 

Wivm= I YimYvmSinOdedct). (CIO) 

JR 
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Then a possible estimator for Ci is given by 



In the limit that Q varies much more slowly than Wi 



U'm 



1 dCl ^ 1 dWll^rr. 



Cl dl Wwm dl' 

we have (7; ~ G- 



(C12) 



1=1' 
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